Elevator-like movements of prestin mediate outer hair cell electromotility.

The outstanding acuity of the mammalian ear relies on cochlear amplification, an active mechanism based on the electromotility (eM) of outer hair cells. eM is a piezoelectric mechanism generated by little-understood, voltage-induced conformational changes of the anion transporter homolog prestin (SLC26A5). We used a combination of molecular dynamics (MD) simulations and biophysical approaches to identify the structural dynamics of prestin that mediate eM. MD simulations showed that prestin samples a vast conformational landscape with expanded (ES) and compact (CS) states beyond previously reported prestin structures. Transition from CS to ES is dominated by the translational-rotational movement of prestin's transport domain, akin to elevator-type substrate translocation by related solute carriers. Reversible transition between CS and ES states was supported experimentally by cysteine accessibility scanning, cysteine cross-linking between transport and scaffold domains, and voltage-clamp fluorometry (VCF). Our data demonstrate that prestin's piezoelectric dynamics recapitulate essential steps of a structurally conserved ion transport cycle.

Apo state pore opening as functional basis of increased EAAT anion channel activity in episodic ataxia 6.

SLC1A2 and SLC1A3 encode the glial glutamate transporters EAAT2 and EAAT1, which are not only the predominant glutamate uptake carriers in our brain, but also function as anion channels. Two homologous mutations, which predict substitutions of prolines in the center of the fifth transmembrane helix by arginine (P289R EAAT2, P290R EAAT1), have been identified in patients with epileptic encephalopathy (SLC1A2) or with episodic ataxia type 6 (SLC1A3). Both mutations have been shown to impair glutamate uptake and to increase anion conduction. The molecular processes that link the disease-causing mutations to two major alterations of glutamate transporter function remain insufficiently understood. The mutated proline is conserved in every EAAT. Since the pathogenic changes mainly affect the anion channel function, we here study the functional consequences of the homologous P312R mutation in the neuronal glutamate transporter EAAT4, a low capacity glutamate transporter with predominant anion channel function. To assess the impact of charge and structure of the inserted amino acid for the observed functional changes, we generated and functionally evaluated not only P312R, but also substitutions of P312 with all other amino acids. However, only exchange of proline by arginine, lysine, histidine and asparagine were functionally tolerated. We compared WT, P312R and P312N EAAT4 using a combination of cellular electrophysiology, fast substrate application and kinetic modelling. We found that WT and mutant EAAT4 anion currents can be described with a 11-state model of the transport cycle, in which several states are connected to branching anion channel states to account for the EAAT anion channel function. Substitutions of P312 modify various transitions describing substrate binding/unbinding, translocation or anion channel opening. Most importantly, P312R generates a new anion conducting state that is accessible in the outward facing apo state and that is the main determinant of the increased anion conduction of EAAT transporters carrying this mutation. Our work provides a quantitative description how a naturally occurring mutation changes glutamate uptake and anion currents in two genetic diseases.

Corrigendum: Diclofenac and other non-steroidal anti-inflammatory drugs (NSAIDs) are competitive antagonists of the human P2X3 receptor.

[This corrects the article DOI: 10.3389/fphar.2023.1120360.].

Domain- and state-specific shape of the electric field tunes voltage sensing in voltage-gated sodium channels.

The ability to sense transmembrane voltage underlies most physiological roles of voltage-gated sodium (Nav) channels. Whereas the key role of their voltage-sensing domains (VSDs) in channel activation is well established, the molecular underpinnings of voltage coupling remain incompletely understood. Voltage-dependent energetics of the activation process can be described in terms of the gating charge that is defined by coupling of charged residues to the external electric field. The shape of the electric field within VSDs is therefore crucial for the activation of voltage-gated ion channels. Here, we employed molecular dynamics simulations of cardiac Nav1.5 and bacterial NavAb, together with our recently developed tool g_elpot, to gain insights into the voltage-sensing mechanisms of Nav channels via high-resolution quantification of VSD electrostatics. In contrast to earlier low-resolution studies, we found that the electric field within VSDs of Nav channels has a complex isoform- and domain-specific shape, which prominently depends on the activation state of a VSD. Different VSDs vary not only in the length of the region where the electric field is focused but also differ in their overall electrostatics, with possible implications in the diverse ion selectivity of their gating pores. Due to state-dependent field reshaping, not only translocated basic but also relatively immobile acidic residues contribute significantly to the gating charge. In the case of NavAb, we found that the transition between structurally resolved activated and resting states results in a gating charge of 8e, which is noticeably lower than experimental estimates. Based on the analysis of VSD electrostatics in the two activation states, we propose that the VSD likely adopts a deeper resting state upon hyperpolarization. In conclusion, our results provide an atomic-level description of the gating charge, demonstrate diversity in VSD electrostatics, and reveal the importance of electric-field reshaping for voltage sensing in Nav channels.

Diclofenac and other non-steroidal anti-inflammatory drugs (NSAIDs) are competitive antagonists of the human P2X3 receptor.

Introduction: The P2X3 receptor (P2X3R), an ATP-gated non-selective cation channel of the P2X receptor family, is expressed in sensory neurons and involved in nociception. P2X3R inhibition was shown to reduce chronic and neuropathic pain. In a previous screening of 2000 approved drugs, natural products, and bioactive substances, various non-steroidal anti-inflammatory drugs (NSAIDs) were found to inhibit P2X3R-mediated currents. Methods: To investigate whether the inhibition of P2X receptors contributes to the analgesic effect of NSAIDs, we characterized the potency and selectivity of various NSAIDs at P2X3R and other P2XR subtypes using two-electrode voltage clamp electrophysiology. Results: We identified diclofenac as a hP2X3R and hP2X2/3R antagonist with micromolar potency (with IC50 values of 138.2 and 76.7 µM, respectively). A weaker inhibition of hP2X1R, hP2X4R, and hP2X7R by diclofenac was determined. Flufenamic acid (FFA) inhibited hP2X3R, rP2X3R, and hP2X7R (IC50 values of 221 µM, 264.1 µM, and ∼900 µM, respectively), calling into question its use as a non-selective ion channel blocker, when P2XR-mediated currents are under study. Inhibition of hP2X3R or hP2X2/3R by diclofenac could be overcome by prolonged ATP application or increasing concentrations of the agonist α,ß-meATP, respectively, indicating competition of diclofenac and the agonists. Molecular dynamics simulation showed that diclofenac largely overlaps with ATP bound to the open state of the hP2X3R. Our results suggest a competitive antagonism through which diclofenac, by interacting with residues of the ATP-binding site, left flipper, and dorsal fin domains, inhibits the gating of P2X3R by conformational fixation of the left flipper and dorsal fin domains. In summary, we demonstrate the inhibition of the human P2X3 receptor by various NSAIDs. Diclofenac proved to be the most effective antagonist with a strong inhibition of hP2X3R and hP2X2/3R and a weaker inhibition of hP2X1R, hP2X4R, and hP2X7R. Discussion: Considering their involvement in nociception, inhibition of hP2X3R and hP2X2/3R by micromolar concentrations of diclofenac, which are rarely reached in the therapeutic range, may play a minor role in analgesia compared to the high-potency cyclooxygenase inhibition but may explain the known side effect of taste disturbances caused by diclofenac.

A Novel Homozygous KLHL3 Mutation as a Cause of Autosomal Recessive Pseudohypoaldosteronism Type II Diagnosed Late in Life.

INTRODUCTION: Pseudohypoaldosteronism type II (PHA II) is a Mendelian disorder, featuring hyperkalemic acidosis and low plasma renin levels, typically associated with hypertension. Mutations in WNK1, WNK4, CUL3, and KLHL3 cause PHA II, with dominant mutations in WNK1, WNK4, and CUL3 and either dominant or recessive mutations in KLHL3. Fourteen families with recessive KLHL3 mutations have been reported, with diagnosis at the age of 3 months to 56 years, typically in individuals with normal kidney function. METHODS: We performed clinical and genetic investigations in a patient with hyperkalemic hypertension and used molecular dynamics simulations, heterologous expression in COS7 cells, and Western blotting to investigate the effect of a KLHL3 candidate disease mutation on WNK4 protein expression. RESULTS: The patient, a 58-year-old woman from a consanguineous family, showed hypertension, persistent hyperkalemic acidosis associated with severe muscle pain, nephrolithiasis, chronic kidney disease (CKD), and coronary heart disease. Therapy with hydrochlorothiazide corrected hyperkalemia, hypertension, and muscle pain. Genetic analysis revealed a homozygous p.Arg431Trp mutation at a highly conserved KLHL3 position. Simulations suggested reduced stability of the mutant protein, which was confirmed by Western blot. Compared with wild-type KLHL3, cotransfection of p.Arg431Trp KLHL3 led to increased WNK4 protein levels, inferred to cause increased NaCl reabsorption via the thiazide-sensitive carrier and PHA II. CONCLUSIONS: Even in patients presenting late in life and in the presence of CKD, PHA II should be suspected if renin levels are low and hyperkalemic acidosis and hypertension are inadequate for CKD stage, particularly in the presence of a suspicious family history.

Molecular Basis of Coupled Transport and Anion Conduction in Excitatory Amino Acid Transporters.

Glutamate is the major excitatory neurotransmitter in the mammalian central nervous system. After its release from presynaptic nerve terminals, glutamate is quickly removed from the synaptic cleft by excitatory amino acid transporters (EAATs) 1-5, a subfamily of glutamate transporters. The five proteins utilize a complex transport stoichiometry that couples glutamate transport to the symport of three Na+ ions and one H+ in exchange with one K+ to accumulate glutamate against up to 106-fold concentration gradients. They are also anion-selective channels that open and close during transitions along the glutamate transport cycle. EAATs belong to a larger family of secondary-active transporters, the SLC1 family, which also includes purely Na+- or H+-coupled prokaryotic transporters and Na+-dependent neutral amino acid exchangers. In recent years, molecular cloning, heterologous expression, cellular electrophysiology, fluorescence spectroscopy, structural approaches, and molecular simulations have uncovered the molecular mechanisms of coupled transport, substrate selectivity, and anion conduction in EAAT glutamate transporters. Here we review recent findings on EAAT transport mechanisms, with special emphasis on the highly conserved hairpin 2 gate, which has emerged as the central processing unit in many of these functions.

Molecular mechanisms of ion conduction and ion selectivity in TMEM16 lipid scramblases.

TMEM16 lipid scramblases transport lipids and also operate as ion channels with highly variable ion selectivities and various physiological functions. However, their molecular mechanisms of ion conduction and selectivity remain largely unknown. Using computational electrophysiology simulations at atomistic resolution, we identified the main ion-conductive state of TMEM16 lipid scramblases, in which an ion permeation pathway is lined by lipid headgroups that directly interact with permeating ions in a voltage polarity-dependent manner. We found that lipid headgroups modulate the ion-permeability state and regulate ion selectivity to varying degrees in different scramblase isoforms, depending on the amino-acid composition of the pores. Our work has defined the structural basis of ion conduction and selectivity in TMEM16 lipid scramblases and uncovered the mechanisms responsible for the direct effects of membrane lipids on the conduction properties of ion channels.

g_elpot: A Tool for Quantifying Biomolecular Electrostatics from Molecular Dynamics Trajectories.

Electrostatic forces drive a wide variety of biomolecular processes by defining the energetics of the interaction between biomolecules and charged substances. Molecular dynamics (MD) simulations provide trajectories that contain ensembles of structural configurations sampled by biomolecules and their environment. Although this information can be used for high-resolution characterization of biomolecular electrostatics, it has not yet been possible to calculate electrostatic potentials from MD trajectories in a way allowing for quantitative connection to energetics. Here, we present g_elpot, a GROMACS-based tool that utilizes the smooth particle mesh Ewald method to quantify the electrostatics of biomolecules by calculating potential within water molecules that are explicitly present in biomolecular MD simulations. g_elpot can extract the global distribution of the electrostatic potential from MD trajectories and measure its time course in functionally important regions of a biomolecule. To demonstrate that g_elpot can be used to gain biophysical insights into various biomolecular processes, we applied the tool to MD trajectories of the P2X3 receptor, TMEM16 lipid scramblases, the secondary-active transporter GltPh, and DNA complexed with cationic polymers. Our results indicate that g_elpot is well suited for quantifying electrostatics in biomolecular systems to provide a deeper understanding of its role in biomolecular processes.

Uncoupling sodium channel dimers restores the phenotype of a pain-linked Nav 1.7 channel mutation.

BACKGROUND AND PURPOSE: The voltage-gated sodium channel Nav 1.7 is essential for adequate perception of painful stimuli. Mutations in the encoding gene, SCN9A, cause various pain syndromes in humans. The hNav 1.7/A1632E channel mutant causes symptoms of erythromelalgia and paroxysmal extreme pain disorder (PEPD), and its main gating change is a strongly enhanced persistent current. On the basis of recently published 3D structures of voltage-gated sodium channels, we investigated how the inactivation particle binds to the channel, how this mechanism is altered by the hNav 1.7/A1632E mutation, and how dimerization modifies function of the pain-linked mutation. EXPERIMENTAL APPROACH: We applied atomistic molecular simulations to demonstrate the effect of the mutation on channel fast inactivation. Native PAGE was used to demonstrate channel dimerization, and electrophysiological measurements in HEK cells and Xenopus laevis oocytes were used to analyze the links between functional channel dimerization and impairment of fast inactivation by the hNav 1.7/A1632E mutation. KEY RESULTS: Enhanced persistent current through hNav 1.7/A1632E channels was caused by impaired binding of the inactivation particle, which inhibits proper functioning of the recently proposed allosteric fast inactivation mechanism. hNav 1.7 channels form dimers and the disease-associated persistent current through hNav 1.7/A1632E channels depends on their functional dimerization status: Expression of the synthetic peptide difopein, a 14-3-3 inhibitor known to functionally uncouple dimers, decreased hNav 1.7/A1632E channel-induced persistent currents. CONCLUSION AND IMPLICATIONS: Functional uncoupling of mutant hNav 1.7/A1632E channel dimers restored their defective allosteric fast inactivation mechanism. Our findings support the concept of sodium channel dimerization and reveal its potential relevance for human pain syndromes.

Single-nucleotide variants in human CD81 influence hepatitis C virus infection of hepatoma cells.

An estimated number of 71 million people are living with chronic hepatitis C virus (HCV) infection worldwide and 400,000 annual deaths are related to the infection. HCV entry into the hepatocytes is complex and involves several host factors. The tetraspanin human CD81 (hCD81) is one of the four essential entry factors and is composed of one large extracellular loop, one small extracellular loop, four transmembrane domains, one intracellular loop and two intracellular tails. The large extracellular loop interacts with the E2 glycoprotein of HCV. Regions outside the large extracellular loop (backbone) of hCD81 have a critical role in post-binding entry steps and determine susceptibility of hepatocytes to HCV. Here, we investigated the effect of five non-synonymous single-nucleotide variants in the backbone of hCD81 on HCV susceptibility. We generated cell lines that stably express the hCD81 variants and infected the cells using HCV pseudoparticles and cell culture-derived HCV. Our results show that all the tested hCD81 variants support HCV pseudoparticle entry with similar efficiency as wild-type hCD81. In contrast, variants A54V, V211M and M220I are less supportive to cell culture-derived HCV infection. This altered susceptibility is HCV genotype dependent and specifically affected the cell entry step. Our findings identify three hCD81 genetic variants that are impaired in their function as HCV host factors for specific viral genotypes. This study provides additional evidence that genetic host variation contributes to inter-individual differences in HCV infection and outcome.

Unique structure and function of viral rhodopsins.

Recently, two groups of rhodopsin genes were identified in large double-stranded DNA viruses. The structure and function of viral rhodopsins are unknown. We present functional characterization and high-resolution structure of an Organic Lake Phycodnavirus rhodopsin II (OLPVRII) of group 2. It forms a pentamer, with a symmetrical, bottle-like central channel with the narrow vestibule in the cytoplasmic part covered by a ring of 5 arginines, whereas 5 phenylalanines form a hydrophobic barrier in its exit. The proton donor E42 is placed in the helix B. The structure is unique among the known rhodopsins. Structural and functional data and molecular dynamics suggest that OLPVRII might be a light-gated pentameric ion channel analogous to pentameric ligand-gated ion channels, however, future patch clamp experiments should prove this directly. The data shed light on a fundamentally distinct branch of rhodopsins and may contribute to the understanding of virus-host interactions in ecologically important marine protists.

Allosteric gate modulation confers K+ coupling in glutamate transporters.

Excitatory amino acid transporters (EAATs) mediate glial and neuronal glutamate uptake to terminate synaptic transmission and to ensure low resting glutamate concentrations. Effective glutamate uptake is achieved by cotransport with 3 Na+ and 1 H+ , in exchange with 1 K+ . The underlying principles of this complex transport stoichiometry remain poorly understood. We use molecular dynamics simulations and electrophysiological experiments to elucidate how mammalian EAATs harness K+ gradients, unlike their K+ -independent prokaryotic homologues. Glutamate transport is achieved via elevator-like translocation of the transport domain. In EAATs, glutamate-free re-translocation is prevented by an external gate remaining open until K+  binding closes and locks the gate. Prokaryotic GltPh contains the same K+ -binding site, but the gate can close without K+ . Our study provides a comprehensive description of K+ -dependent glutamate transport and reveals a hitherto unknown allosteric coupling mechanism that permits adaptions of the transport stoichiometry without affecting ion or substrate binding.

ß1 subunit stabilises sodium channel Nav1.7 against mechanical stress.

KEY POINTS: The voltage-gated sodium channel Nav1.7 is a key player in neuronal excitability and pain signalling. In addition to voltage sensing, the channel is also modulated by mechanical stress. Using whole-cell patch-clamp experiments, we discovered that the sodium channel subunit ß1 is able to prevent the impact of mechanical stress on Nav1.7. An intramolecular disulfide bond of ß1 was identified to be essential for stabilisation of inactivation, but not activation, against mechanical stress using molecular dynamics simulations, homology modelling and site-directed mutagenesis. Our results highlight the role of segment 6 of domain IV in fast inactivation. We present a candidate mechanism for sodium channel stabilisation against mechanical stress, ensuring reliable channel functionality in living systems. ABSTRACT: Voltage-gated sodium channels are key players in neuronal excitability and pain signalling. Precise gating of these channels is crucial as even small functional alterations can lead to pathological phenotypes such as pain or heart failure. Mechanical stress has been shown to affect sodium channel activation and inactivation. This suggests that stabilising components are necessary to ensure precise channel gating in living organisms. Here, we show that mechanical shear stress affects voltage dependence of activation and fast inactivation of the Nav1.7 channel. Co-expression of the ß1 subunit, however, protects both gating modes of Nav1.7 against mechanical shear stress. Using molecular dynamics simulation, homology modelling and site-directed mutagenesis, we identify an intramolecular disulfide bond of ß1 (Cys21-Cys43) which is partially involved in this process: the ß1-C43A mutant prevents mechanical modulation of voltage dependence of activation, but not of fast inactivation. Our data emphasise the unique role of segment 6 of domain IV for sodium channel fast inactivation and confirm previous reports that the intracellular process of fast inactivation can be modified by interfering with the extracellular end of segment 6 of domain IV. Thus, our data suggest that physiological gating of Nav1.7 may be protected against mechanical stress in a living organism by assembly with the ß1 subunit.

Gating Charge Calculations by Computational Electrophysiology Simulations.

Electrical cell signaling requires adjustment of ion channel, receptor, or transporter function in response to changes in membrane potential. For the majority of such membrane proteins, the molecular details of voltage sensing remain insufficiently understood. Here, we present a molecular dynamics simulation-based method to determine the underlying charge movement across the membrane-the gating charge-by measuring electrical capacitor properties of membrane-embedded proteins. We illustrate the approach by calculating the charge transfer upon membrane insertion of the HIV gp41 fusion peptide, and validate the method on two prototypical voltage-dependent proteins, the Kv1.2 K+ channel and the voltage sensor of the Ciona intestinalis voltage-sensitive phosphatase, against experimental data. We then use the gating charge analysis to study how the T1 domain modifies voltage sensing in Kv1.2 channels and to investigate the voltage dependence of the initial binding of two Na+ ions in Na+-coupled glutamate transporters. Our simulation approach quantifies various mechanisms of voltage sensing, enables direct comparison with experiments, and supports mechanistic interpretation of voltage sensitivity by fractional amino acid contributions.

Membrane potentials regulating GPCRs: insights from experiments and molecular dynamics simulations.

G-protein coupled receptors (GPCRs) form the largest class of membrane proteins in humans and the targets of most present drugs. Membrane potential is one of the defining characteristics of living cells. Recent work has shown that the membrane voltage, and changes thereof, modulates signal transduction and ligand binding in GPCRs. As it may allow differential signalling patterns depending on tissue, cell type, and the excitation status of excitable cells, GPCR voltage sensitivity could have important implications for their pharmacology. This review summarises recent experimental insights on GPCR voltage regulation and the role of molecular dynamics simulations in identifying the structural basis of GPCR voltage-sensing. We discuss the potential significance for drug design on GPCR targets from excitable and non-excitable cells.

Structural Mechanisms of Voltage Sensing in G Protein-Coupled Receptors.

G-protein-coupled receptors (GPCRs) form the largest superfamily of membrane proteins and one-third of all drug targets in humans. A number of recent studies have reported evidence for substantial voltage regulation of GPCRs. However, the structural basis of GPCR voltage sensing has remained enigmatic. Here, we present atomistic simulations on the δ-opioid and M2 muscarinic receptors, which suggest a structural and mechanistic explanation for the observed voltage-induced functional effects. The simulations reveal that the position of an internal Na(+) ion, recently detected to bind to a highly conserved aqueous pocket in receptor crystal structures, strongly responds to voltage changes. The movements give rise to gating charges in excellent agreement with previous experimental recordings. Furthermore, free energy calculations show that these rearrangements of Na(+) can be induced by physiological membrane voltages. Due to its role in receptor function and signal bias, the repositioning of Na(+) has important general implications for signal transduction in GPCRs.

Insights into the function of ion channels by computational electrophysiology simulations.

Ion channels are of universal importance for all cell types and play key roles in cellular physiology and pathology. Increased insight into their functional mechanisms is crucial to enable drug design on this important class of membrane proteins, and to enhance our understanding of some of the fundamental features of cells. This review presents the concepts behind the recently developed simulation protocol Computational Electrophysiology (CompEL), which facilitates the atomistic simulation of ion channels in action. In addition, the review provides guidelines for its application in conjunction with the molecular dynamics software package GROMACS. We first lay out the rationale for designing CompEL as a method that models the driving force for ion permeation through channels the way it is established in cells, i.e., by electrochemical ion gradients across the membrane. This is followed by an outline of its implementation and a description of key settings and parameters helpful to users wishing to set up and conduct such simulations. In recent years, key mechanistic and biophysical insights have been obtained by employing the CompEL protocol to address a wide range of questions on ion channels and permeation. We summarize these recent findings on membrane proteins, which span a spectrum from highly ion-selective, narrow channels to wide diffusion pores. Finally we discuss the future potential of CompEL in light of its limitations and strengths. This article is part of a Special Issue entitled: Membrane Proteins edited by J.C. Gumbart and Sergei Noskov.

Molecular physiology of EAAT anion channels.

Glutamate is the major excitatory neurotransmitter in the mammalian central nervous system. After release from presynaptic nerve terminals, glutamate is quickly removed from the synaptic cleft by a family of five glutamate transporters, the so-called excitatory amino acid transporters (EAAT1-5). EAATs are prototypic members of the growing number of dual-function transport proteins: they are not only glutamate transporters, but also anion channels. Whereas the mechanisms underlying secondary active glutamate transport are well understood at the functional and at the structural level, mechanisms and cellular roles of EAAT anion conduction have remained elusive for many years. Recently, molecular dynamics simulations combined with simulation-guided mutagenesis and experimental analysis identified a novel anion-conducting conformation, which accounts for all experimental data on EAAT anion currents reported so far. We here review recent findings on how EAATs accommodate a transporter and a channel in one single protein.